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The electronic orders in Hubbard models on a Kagome lattice at van Hove filling are of intense 
current interest and debate. We study this issue using the singular-mode functional renormalization 
group theory. We discover a rich variety of electronic instabilities under short range interactions. 
With increasing on-site repulsion [/, the system develops successively ferromagnetism, intra unit- 
cell antiferromagnetism, and charge bond order. With nearest-neighbor Coulomb interaction V 
alone (t/ = 0), the system develops intra-unit-cell charge density wave order for small V ^ s— wave 
superconductivity for moderate V, and the charge density wave order appears again for even larger 
V . With both U and V ^ we also find spin bond order and chiral (1^2 _y2 + idxy superconductivity 
in some particular regimes of the phase diagram. We find that the s-wave superconductivity is 
a result of charge density wave fluctuations and the squared logarithmic divergence in the pairing 
susceptibility. On the other hand, the d-wave superconductivity follows from bond order fluctuations 
that avoid the matrix element effect. The phase diagram is vastly different from that in honeycomb 
lattices because of the geometrical frustration in the Kagome lattice. 

PACS numbers: 71.27.+a, 71.10.-w, 64.60.ae, 75.30.Fv 



I. INTRODUCTION 

The Kagome lattice model has attracted considerable 
attention duo to its high degree of geometrical frustra- 
tion. In the Mott insulating limit, several possible states 
have been proposed as the ground state of the Heisen- 
berg model in this lattice, such as the U{1) algebraic 
spin liquid (SL),^ the valance bond solid, ^ the triplet- 
gapped SL,^ and the singlet-gapped SL with signatures 
of Z2 topological order. ^ On the other hand, several ex- 
otic phases have been proposed for the Kagome Hubbard 
model, such as the ferromagnetism at electron filling 1/3 
(or 5/3) per site,^ the fractional charge at 1/3 filling for 
spinless fermions,^ and the Mott transition in anisotropic 
Kagome lattices. ^'^ 

Of particular interest is the possible phases at the van 
Hove filling (the filling fraction is 2/3 =b 1/6 per site), 
where the Fermi surface (FS) is perfectly nested and has 
saddle points on the edges of the Brillouine zone. These 
properties of the normal state makes it unstable against 
infinitesimal interactions. Similar FS appears in triangle 
and honeycomb lattices and were shown to develop, un- 
der short range repulsive interactions, chiral spin-density- 



wave (SDW) state^ or chiral d^.'^ 



idxy super- 



conducting state. ^^'^^ Both states break time-reversal 
and parity symmetries, and are topologically nontriv- 
ial. Given the similar FS, a simple FS nesting argu- 
ment would predict similar phases in the Kagome model. 
This seems to be the case in a recent variational clus- 
ter perturbation theory (with an additional spin disor- 
dered phase). However, as already realized in^^ and em- 
phasized in^^, the interaction vertex viewed in the band 
basis has a strong momentum dependence (matrix ele- 
ment effect). This is because the character of the Bloch 
state on the FS depends on the position of the momen- 
tum. The matrix element effect weakens the nesting ef- 
fect significantly for a local interaction [/, leading to a 
new phase diagram in a recent analytical renormaliza- 



tion group study. Such an analysis would be exact for 
a featureless fermi surface and infinitesimal interactions, 
but its apphcability to the case of finite interactions to- 
gether with perfect fermi surface nesting with van Hove 
singularity is an interesting issue to be addressed. 

The functional renormalization group (FRG) method 
is a differential perturbation theory with respect to the 
increment of the phase space rather than in the inter- 
action itself. It provides the flow of one-particle ir- 
reducible vertex functions versus the running parame- 
ter that controls the phase space. If implemented ex- 
actly the applicability of FRG is not limited by the 
size of the interaction. In practice, however, the ver- 
tices are truncated up to the four-point vertices under 
the assumption that higher order vertices are irrelevant. 
The FRG is promising to address finite interactions and 
treat particle-particle and particle-hole channels on equal 
footing. The applicability of FRG has been demon- 
strated in the contexts of cuprates and iron based 
superconductors.^^ Recently, a singular- mode functional 
renormalization group (SMFRG) method was developed 
and applied to investigate topological superconductiv- 
ity in correlated electron systems with or near van Hove 
singularities. "'^"'^'"'^^ 

In this paper we perform SMFRG study of the model 
at van Hove filling. We discover a rich variety of elec- 
tronic instabilities under short range interactions. With 
increasing on-site repulsion /7, the system develops suc- 
cessively ferromagnetism (FM), intra unit-cell antiferro- 
magnetism (AFM), and charge bond order (CBO). With 
nearest-neighbor Coulomb interaction V alone (U = 0), 
the system develops intra-unit-cell charge density wave 
(CDW) order for small V ^ s— wave superconductivity 
(sSC) for moderate V ^ and CDW appears again for even 
larger V. With both U and we also find spin bond 
order (SBO) and chiral dj;2-y2 -\- idxy superconductivity 
(dSC). Our results are summarized in the phase diagram 
Fig.9. We find that the sSC is a result of CDW flue- 
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tuations and the squared logarithmic divergence in the 
pairing susceptibihty. On the other hand, the dSC fol- 
lows from bond order fluctuations that avoid the matrix 
element effect. The phase diagram is vastly different from 
that in honeycomb lattices. 

The rest of the paper is arranged as follows. In Sec. II, 
we define the model and illustrate the matrix element ef- 
fect. In Sec. Ill, we introduce the FRG method. In Sec. 
IV, we first discuss the leading instabilities at typical 
points in the parameter space, and conclude by a discus- 
sion of the phase diagram. Finally, Sec.V is a summary 
and perspective of this work. 



II. THE MODEL AND THE MATRIX 
ELEMENT EFFECT 
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The Hubbard model we used for the Kagome lattice is 
given by 

H= ^ (4c,a + h.C.) - /iiVe 

^rii^riii ^V^nifij, (1) 

where t is the hopping integral, (ij) denotes bonds con- 
necting nearest-neighbor sites i and j, a is the spin polar- 
ity, /i is the chemical potential, is the total electron 
number operator, U is the on-site Hubbard interaction 
and V is the Coulomb interaction on nearest-neighbor 
bonds. Fig. 1(a) shows the structure of the Kagome lat- 
tice. The different symbols denote the three sublattices, 
and a and b are the two principle translation vectors. 
Fig. 1(b) shows the band structure of the model along 
high symmetry cuts in the Brillouine zone. The lower 
two bands cross at the Dirac point. The highest band 
is a flat band. The dashed line highlights the van Hove 
singularity. Fig. 1(c) shows the normal state density of 
states. The three sharp peaks arise from the van Hove 
singularities in the lower two bands and the third fiat 
band. Fig. 1(d) shows the FS and the character of the 
Bloch states thereon. The FS appears to be perfectly 
nested. However, the character changes along each seg- 
ment. The end points of each segment are saddle points. 
They have pure but different sublattice characters. The 
characters are mixed within the segment as shown by the 
color scale. 

Consider the local interaction U for the moment. Such 
an interaction causes scattering at any wave vector. The 
nested FS would favor scattering connected by the nest- 
ing vectors and would lead to antiferromagnetism in 
usual case. However, since U preserves sublattice in- 
dices, the character variation along the FS causes sig- 
nificant momentum dependence if U is projected to the 
band basis, a matrix element effect as emphasized in^^. 
This effect hampers the scattering significantly. To have 
a better idea of this effect, we calculate the zero frequency 
bare spin susceptibility Xo(q) for site- local spin densities. 



FIG. 1: (a) Structure of the Kagome lattice. Here a = (1,0) 
and b = (1/2, ^3/2) are the lattice unit vectors. The labels 
1, 2, 3 denote the three sublattices. (b) The tight-binding dis- 
persion along high symmetric directions. The dashed line is 
the Fermi level corresponding to the van Hove filling, (c) 
Normal state density of states, (d) Fermi surface and sublat- 
tice weights in the Bloch states thereon. On the endpoint of 
a fermi surface segment the Bloch state is contributed com- 
pletely by one sublattice as indicated by the numbers. Within 
the segment the character are mixed as a superposition of the 
sublattice characters on the two endpoints. The arrow indi- 
cates one of the nesting vectors, (e) The largest eigenvalue of 
Xo(q) as a function of q. 

where q is the momentum transfer. The susceptibility is 
a matrix function in terms of the sublattice labels {a and 

/?), 

Xfi^i) = -^^G"''(k,iw™)G^«(k + q,iu;^), (2) 

k,m 

where T is the temperature, N the number of unit cells, k 
the lattice momentum, Um the Matsubara frequency, and 
G(k, iuJm) the bare Green's function (in the sublattice ba- 
sis). Fig. 1(e) shows the largest eigenvalue of Xo(q) as a 
function of q (for T = O.OOlt). Instead of isolated peaks 
we see branch cuts of maxima in the momentum space. 
These cuts cross at the origin, where there is in fact a 
logarithmic singularity due to the saddle points. (The 
singularity is smeared by the finite size and finite tem- 
perature in the calculation). It is clear that site- local fer- 
romagnetism rather than antiferromagnetism is the most 
favorable spin order, in contrast to the case in the honey- 
comb lattice. ^^'^^ The lesson we learned from the above 
analysis is that for a multi-sublattice system the matrix 
element effect could weaken the nested scattering and al- 
ter the usual intuition regarding FS nesting. There is, 
however, a caveat in this kind of Stoner analysis, since 
it ignores mode-mode coupling between the particle-hole 
channels, and between particle- hole and particle-particle 
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channels. To treat all channels on equal footing we now 
switch to FRG. 



III. THE SMFRG METHOD 

In the following we apply a particular implemen- 
tation of FRG, i.e., the SMFRG, which appears ad- 
vantageous to treat systems with or near van Hove 
singularities. ^^'^^'^^ In this implementation, a generic 
four-point vertex function ri234, which appears in the in- 
teraction cjc2(— ri234)c3C4, where 1 = (k, a) is a dummy 
label indicating the lattice momentum and sublattice la- 
bel, is decomposed into the pairing (P), the crossing (C), 
and the direct (I)) channels as 

r^^^:_k,-p,p+q ^ E/™(k.«'/^)^™n(q)/n(P'^'M), 

mn 

KiZp,Kp+ci -^'^f^{K(^^^)Cmn{ci)fn{p,iy,P), 

mn 

Kl7,P,P+^M ^ E a, i.)i?„„(q)/„(p, ^^, /3). (3) 

mn 

Here fm is a set of ort honor mal lattice form factors. A 
form factor defines a particular composite boson with def- 
inite collective momentum in the particle-hole or particle- 
particle channel, bearing a definite irreducible represen- 
tation under the point group. The fact that the same 
generic vertex can be decomposed into different channels 
refiects the fact that these channels have mutual overlaps. 
The momentum space form factors are related to the real 
counterparts as, /^(k,a,/3) = ErGm /^(r, Q^, /^je"^^-'^ 
where r belongs to a set of bond vectors connecting sub- 
lattices a and P and assigned to m. In our practice the 
bond vectors are truncated up to those connecting the 
eighth neighbors (or third like-sublattice neighbors). In 
the following we use m = to characterize the 

form factor label m, with / indicating the symmetry of 
the form factor, a one of the two sublattice labels, and 
S a basis bond vector that can generate the set of bond 
vectors under the point group. This is applicable since 
we set the symmetry center at an atomic site so that the 
symmetry group is C2. Under this point group, a and f3 
are invariant. There are only two irreducible representa- 
tions Ag (even) and Au (odd) for C2 • We emphasize that 
even though the real-space range of the form factors is 
truncated the range of composite boson scattering is un- 
limited. This enables us to address the thermodynamic 
limit. 

In the SMFRG, P, C and D are substituted into in- 
dependent sets of one-loop and one-particle irreducible 
FRG Feynman diagrams where they would become po- 
tentially singular. (For example P is substituted into the 
particle-particle diagram.) This leads to the differential 
change 9P, dC and dD with respect to the change of 
the running scale A, which we chose as the infrared cut- 
off of the Matsubara frequency. Since there are overlaps 
among the three channels, the full change is a sum of 



the partial one plus the overlaps. It is in this sense that 
SMFRG takes care of mode-mode coupling and treats all 
channels on equal footing. This enables an initially repul- 
sive pairing channel to become attractive at low energy 
scales, and is thus able to refiect the well-known Kohn- 
Luttinger anomaly. The technical details have been ex- 
posed elsewhere. ^^'^^ 

The effective interaction in the superconducting (SC), 
spin density wave (SDW), and CDW channels are given 
by Vsc -P, Vsdw C, and Vcdw = C - 2D, respec- 
tively. By singular value decomposition, we determine 
the leading instability in each channel, 

^r"(q.) = E ^"</'"(m)C(n), (4) 

where x = sc^sdw^cdw, is the singular value of the 
a—th singular mode, (j)^ and t/j^ are the right and left 
eigenvectors of 14, respectively. We fix the phase of the 
eigenvectors by requiring Re[}2^(j)'^{m)2p^{m)] > so 
that < corresponds to an attractive mode in the 
X— channel. In the pairing channel q^c = addresses 
the Cooper instability. The ordering wave vector in the 
SDW/CDW channel 

Q — ^.sdw/cdw is chosen at which 
Vsdw/cdw{^) has the most attractive eigenvalue. We note 
that such a vector has symmetry-related images, and may 
change during the FRG fiow before settling down to fixed 
values. On the other hand, given the most singular mode, 
an effective field can be defined for the ordered state (or 
the condensed composite boson), 

Hsc = ^sc{m)f^{k, a, /3)c^^^cLk,/3,i + 

m,k 

Hcdw= Yl V^cd^M/m(k,Q^,/3)4+q^^^^^^^Ck,/3,a+h.C., 
m,(j,k 

Hsdw = ^V^sd^(m)/4(k,a,/3)c^q^^^^^^^Ck,/3,i + h.c.(5) 

m,k 

up to global factors. It is understood that the sublattice 
labels a and f3 are determined by m according to our 
construction of form factors. The order parameters are 
encoded in the coefficients in the above field operators. 
Two remarks are in order. First there is a residual SU(2) 
degeneracy in the case of triplet pairing and in the SDW 
order parameters. Second, the order parameters are in 
general nonlocal in real space (unless the contributing 
form factors are all local). 



IV. SMFRG RESULTS 

In this section we provide the SMFRG results for the 
model defined in the previous section. We begin by 
discussing the results at specific points in the parameter 
space ([/, y), and summarize the systematic results on a 
dense grid of (/7, V) by a phase diagram. 
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FIG. 2: (Color online) Results for U = 2t and V = 0. (a) 
FRG flow of (the inverse of) the most singular values S in the 
SC (blue dashed line), SDW( green solid line), and CDW (red 
dashed-dot line) channels, (b) The renormalized interaction 
J^m ^sdw foi" ^ — {^9 J tt, 0) (a = 1, 2, 3) as a function of the 
collective momentum q. The hexagon indicates the Brillouin- 
zone boundary, (c) The order parameter hi (drawn as arrows) 
associated with the dominant SDW singular mode. 
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Ferromagnetic order: For U = 2t and V = 0^ Fig. 2 (a) 
shows the flow of the most negative singular values ( 
denoted as S) in the SC, SDW, and CDW channels. 
Clearly the SDW (green solid line) is the leading insta- 
bility. During the flow q^sdw evolves from qi = (tt, 7^/^/3) 
and settles down at q2 =0. The renormalized inter ac- 
tion E™Kdr(q) for ni = {Ag,a,0) {a = 1,2,3), which 
have dominant value in the leading singular mode, is 
shown in Fig. 2(b). It has a strong peak at momentum 
q = 0. Because the dominant form factor is local the 
ordered spin density is site-local. The effective field 
operator Hgdw according to Eq.(5) can be rewritten 
as Hsdw = Xlicr hiacl^Cia- with the order parameter hi 
shown in Fig. 2(c). This describes a FM order. The 
SC and CDW channel turn out to be sub-leading from 
Fig.2(a). 

Intra-unit-cell antiferromagnetic order: For U = 2.5t 
and y = 0, the flow of the singular values is shown in 
Fig. 3(a). Again the SDW channel is the leading instabil- 
ity. During the flow, c^sdw evolves from qi = (7r,7r/>/3) 
and settles down at q2 = 0, in the same fashion as 
above. Fig. 3(b) shows the interaction J2m^sdwi^) 
m = {Ag^a^O) {a = 1,2,3). It also has a strong peak 
at momentum q = 0. There are in fact two degenerate 
singular modes (apart from the SU(2) degeneracy). 
One of them leads to the order parameter hi shown 
in Fig. 3(c), with the ratio : — 1 : 1 on the three 
sublattices. The other mode lead to a ratio 2 : — 1 : — 1 
(not shown). Both modes are antiferromagnetic within 
the unit cell, but is ferromagnetic from cell to cell. 
Comparing to the FM state, we call such a state the 
AFM state, although the ordering momentum is zero. 
A mean field analysis shows that in the ordered state 
the two degenerate modes are mixed in such a way that 
the spin patten is as shown in Fig. 3(d), with an angle of 
120^ between nearby spins. The SC and CDW channels 
remain to be sub- leading from Fig. 3(a). 



FIG. 3: (Color online) Results for U = 2.5t and V = 0. (a) 
FRG flow of the most singular values in the SC (blue dashed 
line), SDW (green solid line), and CDW (red dashed-dot 
line) channels, (b) The renormalized interaction Vldw 
for m = {Ag,a,0) {a = 1,2,3) as a function of q. (c) The 
order parameter hi associated with one of the two degenerate 
SDW modes, (d) The spin structure in the mean field ordered 
state which combines the two degenerate singular modes. 



Charge bond Order: From Fig. 2 and Fig. 3, we find 
that the CDW channel is enhanced with increasing U. 
This trend continues until the CDW channel becomes 
dominant for U > 2.85t. Fig. 4(a) shows the FRG flow for 
U = 3M and V = 0. During the flow the (icdw evolves 
but settles down at q = (0, 27r/V^) (or its symmetric im- 
ages) in the given view field. The dominant renormalized 
interaction Kd^lq)' ^r m = (A,, 1, l/4x + ^3/4^), 
m = (An,2,l/2f) and m = (A^,3,l/4f - V3/4y), is 
shown in Fig. 4(b), where we see isolated peaks at the 
six nesting vectors (three of which are independent 
and correspond to the three form factors). We find 
that the effective field Hcdw constructed according to 
Eq.(5) for the above singular modes can be rewritten 
as Hcdw = T.{ij)aXij{4a^ja + h-c), and is thus a CBO 
state. The pattern of the order parameter Xij depends 
on the ordering vector Q. For Q = (0,27r/\/3), it is 
shown in Fig. 4 (c). Notice that the field Xij is nonzero 
on parallel lines orthogonal to Q. This is also the case 
for the other ordering momenta related to Q by C^y 
operations. Clearly, the CBO breaks both rotation and 
translation symmetries. The reason that the nesting 
vector is at work here is because the bond-centered 
charge density ^^{cl^Cja + cjcr^^cr) connects different 
sublattices, and can take advantage of the inter-saddle 
scattering connected by the nesting vector. Notice that 
this kind of order is already beyond the mean field 
theory. It is a result of the overlap between the SDW 
and CDW channels as seen from Fig. 4 (a) where the 
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FIG. 4: (Color online) Results for U = 3.5t and V = 0. (a) 
FRG flow of the most singular values S in the SC (blue dashed 
line), SDW (green dashed-dot line), and CDW (red solid 
line) channels, (b) The renormalized interaction J2m KXT 
for m = (A^, 1, l/4x + V3/#), m = (An,2,l/2x) and 
m = (An, 3, l/4x — a/3/4^), as a function of q. The three 
independent peaks correspond to the three form factors, re- 
spectively, (c) The real space structure of the order parameter 
Xij associated with one of the dominant CDW modes with the 
ordering momentum Q = (0, 27r/v^)- The red (blue) color in- 
dicates Xij is positive (negative). 



SDW channel dominates at high energy scales, 
pairing channel is still subdominant here. 



The 



Intra-unit-cell Charge density wave: We now consider 
the effect of the nearest neighbor interaction V. Fig. 5 (a) 
shows the FRG flow for /7 = and V = 0.2bt. It 
is clear that the CDW channel (red solid line) is the 
leading instability. During the flow the cicdw evolves 
from qi = (0,2/\/3)7r to q2 = (0,0.385)7r and finally 
settles down at qs = 0. The dominant renormalized 
interaction Em KXT ^ = (^^,^^,0) (c^ = 1,2,3) 
shown in Fig. 5(b) has a sharp peak at q = 0. There are 
two degenerate singular modes. The effective field Hcdw 
constructed according to Eq.(5) can be rewritten as 

Hcdw = EicT Vi^la^if^ with the order parameter r]i shown 
in Fig. 5(c) and (d) for the two singular modes. This is an 
intra-unit-cell CDW state. It beaks rotational symmetry 
but does not break the translation symmetry. It is 
therefore an analogue of the Pomeranchuk instability 
on square lattices. The SDW and SC channels are 
sub-leading in this case. 

S-wave superconductivity: The FRG fiow for C/ = 
and V = OM is shown in Fig. 6 (a). We find that the 
SC channel is the leading instability. Fig. 6(b) shows the 
renormalized interaction VJ^^ for m = (74^^,1,0). (By 
symmetry, interactions involving form factors centered 
on the other sublattices contribute similarly.) Inspection 
of the eigenfunction (j)sc reveals that it has dominant 
values for Ag form factors involving r = and subdom- 
inant values for Ag form factors involving r connecting 
nearest like-sublattice neighbors. The gap function 
from Hsc constructed according to Eq.(5) projected 
on the fermi surface is shown in Fig. 6(c). Clearly it 
is an sSC gap function. Such a pairing symmetry 
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FIG. 5: (Color online) Results for [/ = and V = 0.25t. 
(a) The FRG flow the most singular values S in the SC 
(blue dashed line), SDW (green dashed-dot line), and CDW 
(red solid line) channels, (b) The renormalized interaction 
Em ^clw for m = (A^, a, 0) (a = 1, 2, 3) as function a func- 
tion of q. (c) and (d) show the order parameter r/i associated 
with the two degenerate CDW singular modes. The length 
of the arrows indicate the amplitude and the direction of the 
arrow indicate the sign of the order parameter. 
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FIG. 6: (Color online) Results for [/ = and V = 0.5t. (a) 
The FRG flow the most singular values S in the SC (blue 
solid line), SDW (green dashed line), CDW (red dashed-dot 
line) channels, (b) The renormalized interaction VJ^'^{q) for 
m = (Ag, 1,0) as a function of q. (c) The momentum space 
gap function on the Fermi surface associated with the SC 
singular mode. 



persists for small U > 0. However the dominant pairing 
amplitude for /7 = is on-site, while the amplitude on 
bonds (connecting nearest like-sublattices) increases and 
eventually dominates with increasing U. Inspection of 
Fig. 6 (a) reveals that such an s-wave pairing follows from 
the overlap with the CDW channel. 

Spin bond order. The FRG flow for U = V = 0.75t is 
shown in Fig. 7(a). Clearly the SDW (green solid line) 
is the leading instability. During the flow q_sdw evolves 
from qi = (0,2/>/3)7r to q2 = and flnally settles 
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FIG. 7: (Color online) Results for U = V = 0.75t. (a) FRG 
flow of the most singular values S in the SC (blue dashed 
line), SDW( green sohd hne), and CDW (red dashed-dot 
line) channels, (b) The renormalized interaction VJdw 
for m = (An, 1, l/4x + y3/#), m = (^^,2, l/2x), and 
m = (^„,3,l/4x - V3/4y) as a function of q. The three 
independent peak momenta corresponds to the three form 
factors, respectively, (c) The real space structure of the or- 
der parameter ^ij associated with one of the SDW singular 
modes with the ordering momentum Q = (0,27r/>/3). The 
blue (red) bond indicates that ^ij is positive (negative). 



down at qa = qi. Fig. 7(b) shows the renormalized 
interaction EmKlw m = (A^,l,l/4f + V3/#), 
m = {Au^ 2, l/2x) and m = (A^, 3, l/4x — >/3/4^), where 
we see isolated peaks at the six nesting vectors (three of 
which are independent and correspond to the three form 
factors). The effective field operator in the real space 

can be written as Hsdw = T^{ij)a iij^i4a^j^ + 4^^^^) 
(apart from the SU(2) degeneracy). The pattern of the 
order parameter (^ij depends on the ordering vector Q. 
For Q = (0,27r/V3), it is shown in Fig.7 (c). As in the 
case of CBO state, the order parameter is nonzero on 
parallel lines orthogonal to Q. This describes an SBO 
state. The SC and CDW channel is sub-leading in this 
case. 

Chiral dx2_y2 + idxy superconductivity: Fig. 8 (a) 
shows the FRG flow for U = 2t dind V = IM. 
Clearly, the SC channel is the leading instability. 
Fig. 8(b) shows the renormalized interaction VJ^'^ for 
m = {Ag^ 1, l/2f + ^/3/2y). Such a form factor shows 
the pairing is on third-neighbor bonds (or nearest 
like-sublattice neighbor bonds). From the singular mode 
(psc we construct the effective pairing operator Hgc^ and 
get the gap function in the momentum space as shown 
in Fig. 8(c). This is clearly a dxy-waye gap function. In 
fact there is another degenerate singular mode which 
gives a 6/2^2 -wave gap function (not shown). Using 
the renormalized pairing interaction we performed mean 
field calculations to find that the ordered state is a 



chiral dx2-y2 



idxy superconducting state, which we 



call the dSC state. The chiral state is fully gapped on 
the fermi surface and thus saves more energy. Fig. 8 (a) 
shows that the SDW and CDW channels are dominant 
at high energy scales. Inspection of the later stage of 
the FRG flow reveals that the singular modes in these 
channels contains dominant CBO and SBO components 
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FIG. 8: (Color online) Results for U = 2t and V = 1.5t. (a) 
FRG flow the most singular values S in the SC (blue solid 
line), SDW (green dashed line), and CDW (red dashed-dot 
line) channels, (b) The renormalized interaction Vs^^(q) for 
m = (Ag, 1, l/2x + y/3/2y) as a function of q. (c) The mo- 
mentum space gap function on the Fermi surface associated 
with one of the two degenerate SC singular modes. 



(on nearest bonds). We shall come back to this point 
later. 

The phase diagram: Apart from the typical results 
discussed above, we have performed systematic SMFRG 
calculations on a dense grid in the ([/, V) plane. The 
results are summarized as a phase diagram shown in 
Fig. 9. The CBO and SBO states have ordering momenta 
at one of the nesting vectors, while the others order at 
zero momentum without breaking translation symmetry. 
However, the CDW and AFM states have intra-unit-cell 
structures. This phase diagram can be understood as 
follows. 

Along the U = axis, the s-wave superconductiv- 
ity appears inbetween the intra-unit-cell CDW states at 
small and large values of V. This is counter-intuitive at 
a first sight since increasing V would favor CDW further. 
However, the numerical result is reasonable for the fol- 
lowing reasons. While the CDW susceptibility behaves as 
ln(iy/A) at the running scale A because of the van Hove 
singularities in the normal state density of states (here W 
is of the order of the bandwidth), the SC susceptibility di- 
verges as ln^(W/A) due to a further Cooper instability.^^ 
Therefore, once the initially repulsive pairing channel be- 
comes slightly attractive via the overlap with the CDW 
channel, the pairing interaction could grow in magnitude 
faster than the CDW interaction, and could eventually 
overwhelm the CDW interaction. This explains the emer- 
gence of the s-wave superconductivity for moderate V. 
However, if V is initially small, the overlap with the SC 
channel is small during the flow. On the other hand, if 
V is large enough, the CDW channel diverges before the 
SC channel takes advantage of the fast growth. These 
considerations are consistent with our results along the 
U = axis. 

In the phase diagram we see that both CBO and SBO 
phases are in proximity to the d-wave SC phase. This 
is a reasonable result since we find that the bond orders 
are on nearest-neighbor bonds, while the d-wave pairing 
are on third-neighbor bonds (or nearest like-sublattice 
neighbor bonds). It is the even order processes involv- 
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ing the bond-density interactions that have overlap with 
the above singlet pairing interaction, which are therefore 
immune to the sign structure in the SBO and CBO in- 
teractions. On the other hand, the on-site repulsion dis- 
favors s-wave pairing. This makes d-wave pairing viable. 
Interestingly by utilizing the bond order fluctuations the 
pairing mechanism avoids the matrix element effect that 
would frustrate site-local spin fluctuations at the nesting 
vector. 

Along the V = axis, our SMFRG result predicts the 
charge bond order for large U. This is indeed a spin 
disordered phase as found in^^, and is beyond the mean 
field theory but consistent with the lack of a well defined 
site-local spin ordered phase. The reason that a large 
U favors a spin disordered state rather than local spin 
moment ordering is twofold. First the matrix element ef- 
fect weakens nested scattering and favors ferromagnetic 
ordering. Second a sufficiently large U makes the nested 
scattering more important as compared to the case of 
small U. This would favor antiferromagnetic ordering. 
The site-local spin ordering is thus frustrated by the com- 
petition of ferromagnetism and antiferromagnetism. The 
compromise is the CBO state, which is an analogue of 
the valence bond solid and reflects the short-range spin 
correlations. The limit of U ^ t is beyond the scope 
of our SMFRG, but enables mapping of the model to a 
doped t — J model. We leave it for further investigations. 

Finally, for /7 ^ t and with increasing V, the successive 
orders are FM, CBO, SBO, dSC, sSC and CDW. This 
sequence is reasonable as follows. The CBO and SBO 
states take advantage of V since by connecting different 
sublattices it avoids the matrix element effect. However, 
a large V favors CDW. In the intermediate region, the 
CBO/SBO fluctuations drive dSC while CDW fluctua- 
tions drives sSC, as discussed above. This explains why 
there is a transition from dSC to sSC with increasing V. 

V. SUMMARY AND PERSPECTIVE 

In summary, we have studied the extended Hubbard 
model on Kagome lattice at van Hove filling using 
the SMFRG method. We discovered a variety of 
phases in the parameter space. Along the V = axis 
and with increasing on-site repulsion [/, the system 
develops successively ferromagnetism, intra unit-cell 
antiferromagnetism, and charge bond order. With 
nearest-neighbor Coulomb interaction V alone {U = 0), 
the system develops intra-unit-cell charge density 
wave order for small 5— wave superconductivity for 
moderate V, and CDW appears again for even larger 
V. With both U and V, we also find spin bond order 
and chiral dx2_y2 + idxy superconductivity. We find 
that the 5-wave superconductivity is a result of CDW 
fluctuations and the squared logarithmic divergence in 
the pairing susceptibility. On the other hand, the d-wave 



superconductivity follows from bond order fluctuations 
that avoid the matrix element effect. We summarized 



> 1 



CDW 






sSC 










CBO 








CDW \^^^ 


FM 


afm\^ 



12 3 4 

U 



FIG. 9: The phase diagram of the Kagome lattice at van 
Hove filling. The electronic orders and the associated ordering 
momenta are: FM (Q = 0), intra-unit-cell AFM (Q = 0), 
CBO (Q = M), intra-unit-cell CDW (Q = 0), SBO (Q = M), 
d^i_yi + idxy-v^d^ie SC (dSC, Q = 0) and s-wave SC (sSC, 
Q = 0). Here M is one of the nesting vectors connecting the 
saddle points on the fermi surface. 



the results by the phase diagram in Fig. 9. It is vastly 
different from that in honeycomb lattices, and the 
difference can be attributed to the frustrating matrix 
element effect. 

We notice that the spin 1/2 Kagome lattice has been 
realized in Herbertsmithite ZnCu3{OH)sCl2^^''^^ 
and its isostructural Mg-based paracatamite 
MgCu3{OH)GCl2.^^ Also, the optical Ka gome lat- 
tice has been simulated experimentally in ultra-cold 
atomic gases, and the optical wavelengths can be suitably 
adjusted for fermionic isotopes such as ^Li and ^^i^.^"^ 
With the possibility of tuning U and V continuously, 
the optical lattice with ultracold atomic gases are most 
promising to realize the predictions presented in this 
paper. 

During the writing of this paper we became aware of 
a parallel work in which a similar scenario is addressed. 
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